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ABSTRACT 

The Wavenumber Technique (WT) is a relatively new method 
of underwater sound transmission analysis. One aspect, 
source depth determination, is studied to evaluate its 
validity and test environmental and acoustic sensitivity. 
The horizontal wavenumber spectrum is analyzed to determine 
null spacings in wavenumber space, which indicates source 
depth by the Lloyd’s Mirror interference effect. Comparison 
of this theory with cases of an isospeed sound profile, 
fully absorbing bottom, and flat totally-reflecting surface 
shows excellent agreement for several parametric variations. 
Cases with realistic sound speed profiles and partially 
absorbing bottoms generally agree with theory, but a dis- 
tinct bias is observed. Source depth determination curves, 
which relate the scaled wavenumber spectral intensity null 
spacing to the source depth, are presented for comparison 
with theory and sensitivity analysis. An example is given 
for suggested application of source depth determination. 
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IN TRO DUCTION 



The Wavenumber Technique (WT) is a new method for the 
solution of several problems related to the analysis of 
underwater sound transmission (L auer , 1979) . Its basis is the 
analysis of acoustic propagation in the wavenumber domain, 
which has been described by DiNapoli (1971) and was an 
intermediate step in the computation of transmission loss 
versus range in the Fast Field Program (FFP) . DiNapoli 
(1977) later used the wavenumber domain to study the impe- 
dance of the ocean-bottom interface. Lauer (1979) first 
introduced the WT and its applications to passive localiza- 
tion and tracking and multipath decomposition. He developed 
a range prediction curve based upon peak wavenumbers in the 
transform spectrum for successive range increments from an 
entire pressure range field. 

These investigations presented highly beneficial and 
promising results from the wavenumber depiction, which indi- 
cates that additional research is appropriate. One particu- 
lar aspect of the WT - determination of source depth - will 
be investigated here in order to evaluate its validity in 
comparison with the Lloyd's Mirror effect. Its relation- 
ships to oceanic and acoustic variations will also be 
exa mined . 
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The wavenumber technique as first described by Lauer 
(1979) was based upon input from the FFP, but it has been 
adapted for use with the Parabolic Equation (PE) method of 
propagation loss determination. The PE-based WT is currently 
beinq investiqated by Lauer and others at the Naval Ocean 
Research and Development Activity (NORDA) and was the model 
used in this research at the Naval Postgraduate School 
(NPS) . Descriptions of the PE and WT will be presented to 

display their sensitivities to variations of several model 
parameters. Also, expected relationships of the WT to other 
current propagation prediction models will be discussed to 
consider how results might vary based upon different inputs. 
Particular attention will be given to algorithm differences, 
degrees of approximation, and treatment of the ocean bottom. 

The PS-based WT will then be exercised for several 
parameter variations to evaluate qualitatively the consis- 
tency of the technique. This study will also demonstrate 
the response of the PS-based WT to geometric and actual 
oceanic conditions. Prediction curves for source depth 
determination based upon these variations will be presented 
ana discussed. 
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II. 



WAVENUMBER TECHNIQUE 
A. GENERAL DESCRIPTION 

The WT is a natural alternative for signal depiction 
because it provides information on the directions of energy 
arrivals. It allows reasonably direct physical interpreta- 
tion of ocean acoustic processes and their relationships to 
environmental conditions (Lauer, 1979) . Lauer further states 
that the wavenumber spectrum plot is clearer and more 
informative than typical curves of propagation loss versus 
ranee. The wavenumber K relates the angular frequency and 
the sound speed by 



K = a) = 2irf 

— 3 

c c 



and it has horizontal and vertical components 



K 



+ k z - 



The reference wavenu irber is determined by the sound speed 
minimum 
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The WT will first be described without reference to any 
particular propagation model to present the general nature 
of the algorithm. Fig. 1 is the flow diagram for the WT. 
The left and right sides of this flow diagram correspond to 
the intensity and wavenumber axes, respectively, of the 
spectrum plot that will result from the WT. This figure and 
the following discussion were developed from Lauer (1979 and 
1982a) and augmented by several telephone conversations with 
R. Lauer, T. Lawrence, and R. Evans of NORDA between July 
and December, 1982. 

The WT calculates spectral intensity from the real and 
imaginary parts of .the acoustic pressure 

? = ip(R,Z) e la)t . 



The complex psi versus range field at a specified depth is 
the propagation model product which determines the pressure. 
The complex pressure field can be modified for effects such 
as attenuation, depending upon whether idealized or realis- 
- ic conditions are being investigated. A Fast Fourier 
Transform (FFT) is then applied to yield transformed complex 
pressure. The spectral intensity is de*ermined at each 
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WAVENUMBER TECHNIQUE FLOW DIAGRAM 




Figure 



Wavenumber Technique Flow Diagram. 
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wavenumber increment 



from the transformed complex pressure 



field by 



I * Re T 2 + Im 2 
s T T 



This intensity field is scaled arbitrarily for satisfactory 
display of the relative maxima and minima to define the 
intensity nulls. 

Additional output parameters from the PE which are 
required for the WT include the range step, the number of 
range points, frequency, minimum sound speed and field 
depth. These are the selected variables noted in figure 1. 
The range step is selected to provide sufficient resolution 
to determine the wavenumber null spacing on the spectrum 
plots. The initial radial component of the wavenumber for 
the spectrum is calculated from DiNapoli (1971) by 
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As will be seen in section II-B, if a "scaled” wavenum- 
ber beta is used for the horizontal axis, analysis is sim- 
plified since the nulls in the spectra should be more evenly 
spaced (Lauer, 1979) . The physical meaning of beta will be 
described in the next section. The resulting scaled wave- 
number spectra are then plotted. The wavenumber null spacing 
for each incremental source depth is determined from each 
scaled spectrum. The source depth determination curve is 
generated by plotting these null spacings as a function of 
source depth, as suggested by Lauer (1979). The null spac- 
ing cf an analyzed received signal could be compared with 
this determination, curve to infer the source depth of the 
signal. This application will be discussed later. The 
source depth determination curve will be specific for the 
acoustic and oceanic descriptions of the medium. 

E. PHYSICAL DESCRIPTION 

The preceding description was intended to provide the 
overall concept of the NT whereas whan follows will explain 
rhe physical basis for the technique. The underlying princi- 
ple cf the NT as described by Lauer (1932a) is the Lloyd's 
Mirror effect, which gives the acoustic field by 

p = E(R) e" la3t , 
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where 



E(R) = ex? {iK /r 2 + (Z R - Z $ ) 2 } - exp { iK /r 2 + (Z R + Z g ) 2 } 

/ R 2 + (Z R - Zg) 2 / R 2 + (Z R + Z s ) 2 

for a point source in a semi- infinite medium of constant 
sound speed with a pressure release surface. This is the 
Lloyd's Mirror field in range space, whereas the Lloyd's 
Mirror field in wavenumber space, applicable to this 
research, is given by Lauer (19 82) as 

F(K) = sinCSZg) exp(iBZ^) } 

e 



where beca is defined by 

S = (K 2 -k 2 ) 1 ' 5 

o r > 



\ 



ar.d E (P) and F(K) are related by the Bessel transform pair 



E(R) = /* 2 F(K) J (KR) K dK 
o o 



an d 



-(’<) = 1 /°°E ( R) J (KR) R dR 
2 0 



2h 



For this research the Bessel function was approximated as a 
trigonometric function so that the FFT could be utilized. 
This is considered acceptable since the source and receiver 
are greatly separated, according to Lauer (1982b). Inspec- 
tion of p(K) reveals that the nulls of the spectrum will be 
equally spaced at intervals of 

AS = £ * 

Z S 

This equal spacing in beta is the key element of the appli- 
cation of the WT to source depth determination (Lauer , 1 979) . 

Only the direct and surf ace -ref lectei waves will inter- 
act at a receiver location in a simple illustration of the 
Lloyd's Mirror effect with a smooth surface. This interac- 
tion will be destructive or constructive, as a result of 
their relative phases (Kinsler at al , 1982) and will produce 
nulls or peaks, respectively, in the wavenumber spectrum. 
Fig. 2 depicts the Lloyd's Mirror effect geometry which will 
be used to study test case outputs for null spacing and sub- 
sequent source depth determination. The surface-reflected 
(SR) wave travel distance will increase at a faster rate 
than the direct (D) wave distance as the source is moved 
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Lloyd's Mirror Effect Geomerry. 
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deeper in the water column (i) while the receiver depth (h) 
is stationary. This leads to different arrival patterns and 
nulls at more wavenumbers. Thus the general pattern of 
decreased spacing in beta with increasing source depth is 
observed . 

Bottom interacting waves will be neglected for idealized 
cases. This will be accomplished by the use of a fully 
absorbing bottom. These waves will be important, however, 
in shallow water or at greater ranges. A few test cases will 
have realistic bottoms. Various propagation models handle 
the bottom in different ways, so the inclusion of a bottom 
will be discussed later. Cases in shallow water (lass than 
5000 feet) or an extended ranges (greater than 100 miles) 
will non be considered here. 

Lauer (1 979) observed that the wavenumber spectrum shows 
a series of modulating envelopes which enclose spiices that 
are the eigenvalues associated with the normal modes. Com- 
parison of the wavenumber spectrum with typical transmission 
loss curves reveals that the wavenumber depiction is more 
coherently structured. This leads to a better physical 
understanding of signal transmission because of the rela- 
tionship between horizontal wavenumber and ray theory. He 
gives this as 

kp = cos 9 = 2 Trf cose . 

C 
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The horizontal wavenumber decreases as theta increases. Thus 
the wavenumber axis of the spectrum plot will include wave- 
numbers associated with smaller elevation (or depression) 
angles in increasing order of wavenumber (Lauer , 1 979) . It 
is then expected, for the case of a neglected bottom, that 
spectrum amplitude peaks could be observed that would corre- 
spond to the direct and surface reflected waves. 

C. WAVENUMBER TECHNIQUE APPLICATION 

Lauer (19 7 9) recommends a receiver that is a single 
omnidirectional hydrophone and a source-generated continuous 
wave (CW) received signal for source depth determination. 
He further states that in actual use the single hydrophone 
could be replaced by an array to improve the signal-to-noise 
ratio and to obtain bearing and bearing rate information. 

An example of a possible application of the WT, source 
depth determination, is presented in Fig. 3. A propagation 
model could be run for a series of incremental source depths 
to yield the null spacing relationship, as was shown in Fig. 
1. That curve would be specific to the properties of the 
ocean and values for parameters such as frequency and 
receiver depth. A received signal as shown on the right side 
of Fig. 3 could then be analyzed. Its null spacing could be 
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UT APPLICATION: SOURCE DEPTH DETERMINATION 




igure 3. 



WT Application for So urea Depth Determination. 
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compare! to the model-produced curve to infer the depth of 
the source. 

D. TECHNIQUE LIMITATIONS 

The success of the WT is related directly to the accu- 
racy of the predicted pressure field produced by the propa- 
gation model. There are several possible differences 
between propagation models which might affect the pressure 
field. Some of the more significant include approximations 
of the wave equation, bottom effects, initial source func- 
tions and input limitations. These will be compared for sev- 
eral models to determine relevant restrictions. The pressure 
field will also respond to changes in the ocean. A series 
of different scenarios will be run to provide an initial 
estimate of the sensitivities to some of these factors. 

Any additional limitations will be related to computer 
processing rime and storage requirements. These factors will 
not be addressed directly in rhis research because the pro- 
cessing following the propagation model is minimal. There is 
also no requirement rc store the complex pressure fields or 
rheir transforms once the null spacing has been determined. 
Thus the WT computer utilization is directly proportional to 
the time and storage required for the propagation models. 
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PROPAGAT ION MOD ELS 



A. INTRODUCTION 

Modeling the propagation of sound in the ocean is com- 
plicated by the variability of the ocean, the great range of 
frequencies of interest, and the many applications of sound 
propagation (DiNapoli and Deaven port , 1979) . It is not sur- 

prising that no single model is presently capable of 
addressing all of these variations; there are many semi-re- 
strictive methods. In general all models consider the sound 
speed to be a function of the spatial coordinates and inde- 
pendent of time (DiNapoli and D eavenport, 1979) . Only the 
frequency, water column geometry and oceanic description 
will affect propagation. 

All models, regardless of application, can be divided 
into two classes: range independent or range dependent. 

Range independent models assume that the ocean is cylindri- 
cally symmetrical, the speed of sound is a function only of 
depth and all boundaries are planar perpendicular to the 
depth axis (DiNapoli and Deavenport, 1979) . Range dependent 
models can provide better approximations of real conditions 
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by allowing the sound speed to be a function of two or three 
spatial variables. They also do not require planar bounda- 
ries perpendicular tc the depth axis, thus cylindrically 
symmetric bottom topography can be included (DiNapoli and 
Deavenport , 1 979) . A range-dependent model, the split-step 
Fast Fourier Transform PE (SSFFT) was the numerical model 
used in this research to generate the pressure fields for 
the wavenumber analyses. This PE will be qualitatively com- 
pared to two other ra nge- depends nt algorithms, a range-modi- 
fied Normal Mode (NM) model and a Finite Difference (FD) 
model. It will alsc be compared tc a range-independent 
method, the Fast Field Program (FFP) model. The range-de- 
pendent capability cf the SSFFT was not used in this 
research since this was an initial invest igation. Its com- 
parison with both range dependent and independent models 
may, however, indicate possible model effects on the WT 
which could be studied later. 

In general ail models attempt to solve the acoustic wave 
equation. This equation is valid for sound propagation in 
fluid media with specified boundary conditions 
(Dr ick , 1 97 5, p . 1 1 4) . It becomes the reduced elliptic Helm- 
holtz equation in spatial coordinates under the assumption 
of harmonic time dependence. 
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B. PARABOLIC EQUATION 
1 • Intro d uc ti on 

The PE algorithm used in this research was developed 
by 3rpck (1 978). The basic feature of this PE is the 
replacement of the reduced elliptic Helmholtz equation by a 
parabolic partial differential equation that utilizes the 

Tappert- Hardin split-step Fourier algorithm (SSFFT) for the 

/ 

numerical integration (Brock, 1978) . The original computer 
code as described by Brock has been modified to incorporate 
the source and attenuation descriptions ot Tatro (1977) and 
to include an ocean bottom by Stieglitz eh al (1979). This 
is the version of the PE that is currently available at NPS 
and was the code (PEMODEL) used for this research. (It must 
be emphasized tha 4 - some difficulty was encountered in trying 
no ascertain the correct description of rhe NPS version 
because of inadequate and erroneous comment card documenta- 
tion in the source cede and incomplete manuscript descrip- 
tions.) 

This section will discuss the SSFFT, the original 
Brock algorithm, modifications to that code, and subsequent 
evaluations of the SSFFT. Particular attention will be 
given to these aspects which are expected tc affect the NT. 
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2. P arabo lic App ro x imat ion 



The parabolic approximation to the wave equation was 
first introduced by Leontovich and Fock (1946) as a solution 
for the propagation of electromagnetic waves along the sur- 
face of the earth. Its first introduction to underwater 
acoustics was by Hardin and Tappert (1973) who used the 
split-step technique in conjunction with the FFT to solve 
this particular approximation to the wave equation. The 
result was the SSFFT. Tappert (1977) related the use of the 
parabolic approximation to sound channel propagation in a 
waveguide that is thin vertically and elongated horizontally 
to the range of the first convergence zone or farther. 

Long-range propagation, with which the SSFFT is 
basically concerned, will usually be for low frequencies 
because volume absorption increases strongly with frequency 
(Tappert , 1 977) . The maximum elevation (and depression) 
angles of propagation for long range transmission must be 
small to satisfy the PE (Tapper t, 1977) . The derivation of 
the SSFFT will be discussed later in a simplified version as 
presented by Brock (1 978). 

It appears that the SSFFT is intended for use in 
situations where bottom interaction and surface scattering 
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would be negligible. Thus the use of the SSFFT to examine 
the Lloyd's Mirror effect in the real world may be inappro- 
priate if the ocean surface is sufficiently rough; addi- 
tional errors may occur in situations where bottom effects 
are significant. 

3. Broc k A lgo rithm 
a. Description 

The SSFFT is valid for acoustic pressure in a 
medium of constant density with a monofrequency source and 
cylindrical symmetry about the depth axis (Brock , 1 978) . The 
acoustic pressure can be written as 

p(R,Z) = i|»(R,Z) H 1 ( K R) 

where the zero-order Hankel function of the first kind 
relates the acoustic pressure to ar. outward propagating 
cylindrical wave envelope function (Brock, 1 978) . This is 
valid because of the assumption that at low frequencies all 
significant energy will propagate approximately horizontally 
away from the source (Tappert , 1 9 77) . The asymptotic form of 
the Hankel function is 

H 1 ( K R) = (2/ttK R}* 5 exD { i ( K R - tt / 4 ) } , K R >> 1 . 

0 0 0 * 0 0 • 
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if the receiver is many wavelengths from the source. This 
leads to (Brock, 1978) 

^ RR + 2iK Q <j> R + ^ zz + K Q 2 {n 2 (R,Z) - 1}* = 0 . 

The additional assumption of neglecting the far field 
effects , 




< < 




5 



is the "parabolic" approximation for radial transmission 
(Brock, 1 97 8) . The result is the Leontovich-Fock (1946) PE 



= i ( A+3 ) i|) 

I\ ^ 



where 



A. = -_1 3_^_ 

2K 3z 2 
0 



and 



3 = Ko^ (n 2 - 1) 
2 
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which is the form used by Brock (1978). 
b. Implementation 

Brock (1 978) used the SSFFT to solve the PE 
because it has several significant advantages which outweigh 
its disadvantages. Advantages include exponential accuracy 
in depth, second-order accuracy in range, energy conserva- 
tion, unconditional stability and computational efficiency. 
Disadvantages are a uniform mesh and periodic boundary con- 
ditions to satisfy the FFT, and filtering of the sound speed 
profile by smoothing discontinuities to avoid spurious high 
angular- freguency components (Br ock, 1 979) . Additionally the 
algorithm assumes a flat pressure-release surface, a vanish- 
ing field at the maximum depth and a pseudo-radiation condi- 
tion at the water-bottom interface by smoothly attenuating 
the field (Brock , 1 978 ) . 



A numerical algorithm cf 



= i(A+3) 'b 

i\ 



if/(R+AR,Z) = e 



iAR( B+A) 



«(R,Z) , 



leads to 



y(_R+AR,Z) - e 



iARB 



FFT 1 { e 



iARA 



FFT (. (. R , Z ) ) } 
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or 



iARK ( n 2 -l ) / 2 iARK 2 /2K 

ip(R+AR,Z) = e 0 FFT” 1 {e °FFT( i|»(R ,Z) ) } 



to facilitate the computations since a spatial FFT and its 



inverse are required to transform between depth and wavenum- 



ber (Br oc k, 1 978) . 


This 


procedure is implemented 


by 


t wo 


alternating steps. 


t be 


first of which considers propagat 


ion 


in a homogeneous 


me d ium 


to account for diffraction 


and 


the 


second to account 


f or 


refraction (Brock, 1 978) . 


In 


his 


development Brock 


gi ves 


an alternate calculation 


of 


the 



envelope function stemming from 



1>(R+AR,Z) 



iARA/ 2 iARB iARA/2 
e e e 



R , Z ) t 



This differs slightly from the previous expression contain- 
ing the FFT and inverse FFT. Thus it is not currently known 
which relationship is used by Brock; this could be resolved 
later by a detailed examination of the computer code. In 
any event as McDaniel (1975a) has shown, either approach 
would be of sufficient accuracy with appropriate selection 
of AR and AZ. 
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The SSFFT appears ho be guite satisfactory 
numerically, but it has to simplify the actual oceanic vari- 
ability to allow realistic computation. It apparently will 
handle source and receiver depths and horizontal separation 
distances well as noted by its degrees of accuracy in depth 
and range. The smoothing of the sound speed gradients and 
bottom effects will have to be considered. 

The intent of this description has been to pro- 
vide a base from which comparisons can be made with other 
numerical solutions and modifications to this algorithm. 
The effect of the source model on the WT can be evaluated by 
these comparisons. The primary elements of the WT are the 
real and imaoinary parts of the complex pressure amplitude. 
Thus the algorithmic differences with respect to computation 
of the pressure field will be the primary basis for the com- 
parison of effects on the WT. 

McDaniel (1975a) compared the basic Tappert ana 
Hardin SSFFT, which is used by 3rocx, to normal mode theory. 
She found -hat errors arise in phase and group speeds, 
although one mode can be propagated with the correct ampli- 
tude and group and phase speeds. McDaniel describes three 
sources of error: (1) approximating the wave equation by the 
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SSFFT as no- ‘rd above, (2) limiting the range steps and per- 
missible sound speed gradients, and (3) truncating the field 
at a finite depth. The result of the SSFFT will be a shift 
in the modal interference at long ranges if many modes are 
propagating (McDaniel , 1975a) . The error due to range step 
limitations is third order in the range step increment 
(McDaniel, 1975a) , thus a small range step is required to 
minimize this error. The third error is related to inclusion 
of a highly absorbing bottom which causes preferential 
attenuation of higher-order modes; the dominant contribution 
will thus result from the lower modes (McDan iel, 1975a) . 

Brock et al (1977) have described a procedure 
to improve this condition. It maps the index of refraction 
and depth such that the phase speeds and turning points of 
the normal modes are preserved. This involves construction 
of a "pseudoproblem" such that the parabolic phase speeds 
will be equal to the elliptic phase speeds of the corre- 
sponding modes in the original problem (Brock et al , 1977) . 
This is giver, as 



( n ,Z) 



(n,Z) 



{ ( 2n-l 
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It is found to greatly improve the SSFFT agreement with the 
elliptic solution. This correction technique has one basic 
constraint, which is that isospeed regions are not allowed 
in the sound speed profile because of poor eigenfunction 
mapping (Brock et al, 1977) . 

4. Subseq uent Modifications 

a. Introduction 

There have been two significant modifications 
since the original algorithm was published by Brock (1 978). 
These are source function and volume attenuation alterations 
(Tatro, 1 977) and the incorporation of a bottom and variable 
range step (Stieglitz ex al , 197 9) . The Tatro modification 
was developed in 1977 but was added to the 3rock algorithm 
after 1978. Discussion of these two aspects will complete 
the description of the PE algorithm installed at NFS and 
utilized in this study. 

b. Tatro Modifications 

The Gaussian source originally utilized by Brock 
(1978) was considered to be inefficient in terms of program 
size limitation and an upper frequency limit which is rela- 
tively low (Tatro, 1977) . A new source function was devel- 
oped. This was intended to solve the problem of large 
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vertical wavenumbers from a Saussian source that were previ- 
ously eliminated with range (Tatro, 1977) and led to an 
incorrect depiction of the field. 

This source is basically a low pass filter in 
-he vertical wavenumber domain which is centered at the 
source depth (Tatro, 1 977) . The initial vertical wavenumber 
field will thus be a constant u p to a prescribed value and 
zero for higher wavenumbers. The use of the filter is 
claimed to yield improvement over the Saussian source by 
minimizing the aliasing which might occur upon transforming 
to the depth domain (Tatro, 1977) . The benefits derived from 
the use of this new source are reduced computation time and 
storage and the availability of higher frequencies 
(Tatro, 1 977) . 

He also modified the Brock algorithm by intro- 
ducing volume absorption as calculated from the equation of 
Thorp (1 967) . The decision can be made whether or not to 
include volume absorption. 

c. Stieglitz et al Modi f ica tions 

Brock (1 978) had considered the ?E most useful 
for long range propagation at low frequencies and small 
grazing angles. The bottom was modeled as fully absorbing so 
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that the field would vanish at the maximum depth of the 
transform. This depth was obtained by extending the water 
column depth by one- fourth to obtain the transform depth. 
The index of refraction in this region was 



n 



2 _ 



= n. 



lae 






max 



where a and b are empirical constants. The energy in the 
bottom was attenuated such that any additional modes which 
resulted from transform truncation were minimized. This 
empirical expression was selected to coincide with a compa- 
rable normal mode solution. This rather crude model of the 
bottom appeared to be acceptable for situations wherein the 
bottom had a relatively small effect. 

Stieglitz et al (197 9) considered it necessary 
-c attempt to develop a more realistic ocean bottom since 
acoustic propagation was being utilized in more varied and 
extreme situations where bottom interaction was important. 
They developed two additional options for handling the bot- 
tom: (1) a partially absorbing bottom with specification of 
reflective loss versus grazing angle of the equivalent ray 
and (2) the direct insertion of sediment sound speed and 
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attenuation profiles. The algorithm derives the sediment 
profiles from the grazing angle loss values in the first 
option. The supplied profiles are used directly in the sec- 
ond option (Stieglitz et al ,1979). A new complex modified 
index of refraction 



M(R,Z) = { n 2 ( R , Z ) - 1} + i2b(Z)/K 

o 

is computed if the bottom attenuation is to be calculated. 

One additional change, the calculation of the 
range step, has also been added by Stieglitz et al (1979). 
It consists of an initial search of the vertical wavenumber 
domain for the maximum pressure component. This is then fol- 
lowed by a search for the first component which is 50 dB 
below the maximum component. The wavenumber of this 

com pcnen t 



k = K sine 
z o 50 



is stated to require a range step of 

AR = A /(I - cose ) , 
0 5 0 
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where 



A = 2u/K . 

0 0 

The range step was fixed in the current research for the 
benefit of the spatial transform in the WT, thus the above 
description was not utilized. It has been included to com- 
plete the description of the computer code which apparently 
comprises the NPS model, PEMODEL. 

5. SSFFT E valua t ions 

The SSFFT has been investigated on several occasions 
to examine specific items which might be modified to yield 
better agreement with actual propagation. Most of these 
studies deal with the algorithmic description by Tappert and 
Hardin ( 1973). Examination of these studies will note fur- 
ther possible limitations on accuracy in the SSFFT which can 
be extended to the WT . 

Fitzgerald (1 975) examined the accuracy of the SSFFT 
as a function of range and frequency. It was shown that the 
SSFFT demonstrated greater accuracy at longer ranges, the 
lower the frequency. The SSFFT was considered a good 
approximation in the examples given for 100 HZ and 10 HZ to 
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about 111 and 15000 km, respectively (Fitzgerald, 1975). 
The current study will examine transmission at 50 and 100 HZ 
at ranges to 100 nautical miles. These ranges will exceed 
the limits of Fitzgerald in some cases. 

McDaniel (197 5b) investigated splitting the total 
field into transmitted and reflected fields in deriving the 
SSFFT. It was necessary to neglect the reflected field for 
the matrix associated with the Tappert-Hardin SSFFT, since 
the two fields did not decouple where the wavenumber was 
independent of range (McDaniel, 1975b) . Thus an exact 
expression could not be determined for the transmitted field 
without the reflected field. (The reflected field is neg- 
lected in the current research since it uses the SSFFT.) 
This will lead to phase speed and group speed errors as 
McDaniel (1975a) pointed out earlier, as discussed in Sec- 
tion III-B-4. These errors were minimized by Brock et al 
(1 977) by a mapping technique for both the index of refrac- 
tion and the depth. Even so, the implication of this possi- 
ble error cz the WT is that the interference patterns may 
r.ot be exact, which will lead to inexact or biased wavenum- 
ber nulls. McDaniel (1 975b) proposed another matrix which 
achieved the desired decoupling and thus allowed expressions 
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for the transmitted and reflected fields. This alternative 
was not investigated in this initial research. It may be a 
way to resolve non-exact results and could be studied later 
to isolate quantitatively the biases related to various 
input parameters for the WT. 

Palmer (1976) related possible inaccuracies in the 
SS? FT to approximations involving relatively small terms in 
the exponential approximations for the operators A and B and 
also to the use of a single reference wavenumber. These 
approximations will net be investigated specifically in this 
research but they may be related to potential disagreement 
between the SSFFT and the theoretical Lloyd's Mirror effect. 
The association of errors due to a single wavenumber can be 
recognized from the wavenumber spectra, which will be pre- 
sented later in section IV-A. They will show not only an 
energy concentration at the reference wavenumber but signif- 
icant energy distributed, as well, over a range of values. 
Palmer (1976) emphasized that the physical properties of the 
ocean must be considered before terms can be neglected, such 
as the spreading of energy over several wavenumbers. 

DeSantc (1977) studied the relationship between the 
acoustic pressure and the envelope function (called velocity 
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potential by DeSanto), which is directly relevant to the 
current research where pressure is calculated from the mod- 
eled envelope function. DeSanto arrives at 

4>(R,Z) = A q R 1 Za f ? ( t , Z ) Q(R,t,Z) exp ( iK Q /Zt ) ( R 2 + t 2 ) t a ~ 3/ ^ 2 dt 
where Q satisfies 

Q RR + 2{(1 ~ a)/R + (iK Q R)/t}Q R + Q zz + 2iK Q Q t + 2Q z 3_{ln( p( t ,Z) ) } 

3 z 

K q 2 {K 2 ( t,Z) - K 2 (R,Z) }q 

and 

K 2 (R,Z) = K(R,Z) - K 1 (Z) 4t _ 



The PE for the envelope function is obtained if this func- 
tion is evaluated by stationary phase; thus, it appears that 
a more accurate result could be obtained by better evalua- 
tion of the integral. DeSanto et al (1978) developed a cor- 
rection for the SSFFT, known as the Corrected Parabolic 
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Approximation (CPA), that handled phase inaccuracies better 
than the SSFFT. Further study might address the response of 
the KT to the CPA with respect to improved phase accuracy as 
a function of range. 

6. Environmenta l Sensitivity of the PE 

This final section in the analysis of the PE will 
consider the expected environmental sensitivity of the PE, 
and, by extension, the WT. The important relationship of the 
PS to the WT is the input of the complex pressure amplitude 
and ultimately the resultant wavenumber distribution. As 
discussed earlier in Section II- B, the wavenumber could be 
associated to the source angle and, by extension, to the ray 
type or family (Lauer, 1979). Consequently, oceanic parame- 
ters which modify the transmission of certain rays can be 
related to the resulting wavenumber distribution. 

The sound speed profile will favor certain ray paths 
at different depths in the water column. There should be a 
corresponding shift in ray paths and hence wavenumbers as 
the source and receiver depth geometry is changed. The 
direct and surface-reflected rays will travel at different 
speeds depending upon the gradients along the profile and 
thus lead to different spatial interference patterns. This 
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can leai to wavenumber shifts and changes in the null spac- 
ing. Profiles used in the current study will include real 
examples and isospeed cases. 

Another important aspect is the transmission geom- 
etry. Generally the model will be run at considerable range 
and in relatively deep water, which will be consistent with 
the applicability of the PE. A series of source depths with 
a single fixed receiver depth will be run to generate the 
wavenumber null spacing values that will formulate a deter- 
mination curve for each case of geometrical variation. The 
distribution of these resulting curves can be examined to 
evaluate the effects cf transmission geometry on the WT. 

Two frequencies will be examined in a similar man- 
ner. This would be one cf the most important parameters to 
investigate in future studies since it would have a direct 
bearing on operational applications. It is expected that 
there will be differences in the resulting wavenumber spec- 
tra depending on the input frequency. The relationship 
between frequency and transmission geometry will also be 
studied since these are the easiest parameters to vary dur- 
ing an experiment and will lead to the most complete analy- 
sis cf model variability. 



50 



The surface will be assumed to be planar and pres- 
sure-release so that there will be complete specular reflec- 
tion. The minimum source depth used will be that at which a 
definable interference pattern can be observed for the spec- 
ified frequency. The slope approaches zero near the surface 
or. the source depth determination curves. This result may 
indicate that the WT in its current form is probably not 
applicable for sufficiently shallow sources. 

The final important area of variability will be the 
bottom boundary conditions. There are several options now 
available in the PEMODEL. The fully absorbing and partially 
absorbing bottoms will be considered. A fully absorbing 
bottom will be used in most cases to simplify the interfer- 
ence pattern by eliminating any rays which interact with the 
bottom. Thus the basic interference pattern may be evalu- 
ated. Various bottoms can then be tested to study the 
effect of the bottom on the WT. Differing bottom less 
curves will be compared in a few cases to demonstrate quali- 
tatively that the bottom does have an effect. The third bot- 
tom option of a realistic sediment sound speed could also be 
investigated once the general affect of the bottom has been 
determined. It would not offer any additional significant 
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results at this point beyond those realized from the fully- 
and partially- absorbing bottoms. 

All input parameters are related and thus the exami- 
nation of any one particular factor must be considered with 
respect to the full range of variation of the other factors. 
The number of combinations of parameters would be very large 
for a model such as the PEMODEL. Thus this research can 
only begin to explore general relationships. Additional 
research is necessary to delineate specific parameter 
impacts on the WT. 

C. FAST FIELD PROGRAM 

1 . Genera l Descr ip t i on 

The Fast Field Program (FFP) was selected as the 
first model to compare to the PEMODEL since it was the ori- 
gin of the WT and the initial testing of the WT by Lauer 
(1979) was performed with the FFP. The description of the 
FFP which follows is summarized from DiNapoli (1971). 

The basic result of the FFP is direct numerical 
integration by the application of the FFT to field theory to 
compute propagation predictions in a minimum computation 
time. The pressure can be represented in the form 

o(R,Z) = r GCZ-r, ,Z C ;K) H 1 (KR) K dK , 

- 00 R S o > 
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where G r the Green’s function, must satisfy 



d 2 G + { K 2 ( Z ) -k 2 }G = -6 ( Z-Z ) 
r s 

dz 2 



and the associated boundary conditions. The field equation 
can then be written as 



p(Z,R ) 
m 



/p” AK(2/ni)^ e^o^m /R ^ 
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imR AK 
e o 
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where 



TT 



G(Z,Z C -,K ) K ** 
’ S ’ m m 



imR AK 
e o 



In this form, ? can be obtained through use of the FFT at 
each range increment. Given p as a function of range, the ST 
could be applied as usual. DiNapoli compared an example 
from the FFP with the normal mode program of Bartberger and 
Ackler (1973) and found excellent agreement between the 
locations of the peaks and the real parts of rhe eigenval- 
ues. The only significant difference is that the FFP 
includes more higher order modes bur does not include the 
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first four modes. DiNapoli states that the FFP offers a 
significant reduction in computation time yet still provides 
reasonable results as compared to normal mode theory. This 
is accomplished by expressing the sound velocity profile as 
one or more exponential functions of depth. These exponen- 
tial functions allow the use of recurrence relations to 
quickly calculate the input to the FFT. 

The FFP offers two options for the modeling of the 
bottom; a two-layered fluid bottom with specified sound 
speeds and a semi-infinite bottom of one sound speed. The 
propagation effect from the first would be eventual return 
of rays once reflected from the second bottom interface. The 
semi-infinite bottom would not allow any return of sediment 
rays. 

2. Compa rison of FFP and PE 

The FFP is a range independent model and the PE is 
range dependent. The environment was considered constant 
ever the range in the current research, and thus the two 
models can be compared. The FFP uses the sound speed profile 
as an exponential function wnich is a different approach 
from the linear segmentation of the PE. The number of expo- 
nential functions could be sufficient to model an actual 
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profile approximately as well as the PE. The PE bottom mod- 
eling does net appear to offer any significant advantage in 
theory over the FFP, which allows for a semi-infinite or 
layered bottom. It is likely from a theoretical standpoint 
that the PE and FFP inputs to the WT could be similar enough 
to allow comparison. The wavenumber spectrum shows a series 
of peaks which result from the nearby singularities of the 
Green's function. These can be expressed as the normal 
modes. The utility of the WT stems from the spacing of the 
nulls between the envelopes which contain these peaks. 

D. NORMAL MODE THEORY 

1 • G e ne ral Desc r option 

The fundamental importance of the normal mode solu- 
tion, and the reason that is it used to check the validity 
of other models, is that it is an exact solution to the wave 
equation (Kir.sler et al , 1982 ,p. 430) . The following descrip- 
tion continues from K ins le r et al (1982, pp. 430-432) . 

The solution for a point source is given by 

{7 2 -3 2 /C 2 3t 2 ) p = 1 5($-$ o ) e^ ut , 

4 TT 

where 

5 = (0,Z C ) . 

o S 
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Analysis yields 



p = e Jt0t (-jir)2Z(Z c; )Z(Z)H H K R) 
n n b n o n 

and Zn must satisfy 

d 2 Z + {w 2 /C 2 (Z) - K 2 }Z = 0 
n n n * 

dz 2 

where Kn is constant and appropriate boundary conditions. 
This solution is for trapped modes and does not include 
those which are evanescent. Normal mode theory can be 
linked to the WT by the spatial representation of the wave- 
number 

K(Z) = to 

1 ■■ - ' • 

C(Z) 



The rays in the deep sound channel which have the same local 
direction of propagation correspond to the normal modes by 
the relationship 



cos 0 = 

K(Z) 
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which extends to the normal mode phase speed 




where 1 designates the 1th mode. 

Normal mode theory lends itself to many applications 
in underwater acoustics and it may be formulated for both 
range dependent and independent situations. A dependent form 
which is appropriate for comparison with the SSFFT was 
developed by Kanabis (1975). He developed a normal-mode 
model which allows large changes in the sound speed profile, 
depth, and bottom composition with range. The range domain 
is segmented into regions within each of which sound speed 
and bottom composition are functions of depth but not range. 
The normal mode solutions for each region are matched at the 
interfaces . 

2. Co mpar ison of NH a nd PE 

It would be expected that the WT would respond simi- 
larly to these two models. This is based upon this example 
of a normal mode program and supported by the previously 
discussed evaluation studies involving PE and NM. The NM 
solutions are generally applicable to long range propagation 
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at lew frequencies, which is quite similar to the PE. The 
NM can be identified closely with the WT in a physical 
sense, because of the interference effects in the NM that 
result from several multiple reflections (Officer, 1958, p. 
117). Officer further states that the number of reflections 
will increase and the time intervals between the incidence 
of successive reflections will decrease as the range 
increases. This might be comparable to depth variations 
which will be shown for the WT. The NM yields exact solu- 
tions for modeling with planar boundaries. Thus the results 
of a WT based on NM would be limited only by approximations 
which were intrinsic to the WT. A comparison of PE- and NM- 
based WT data should aid in evaluating the impact of the 
particular PE algorithm. 

E. FINITE DIFFERENCE METHODS 
1 • General Desc r iptions 

One final group of methods which has evolved 
recently are the finite difference methods of solution for 
the parabolic equations (Lee and Papadakis, 1979) . They 
seek a more general solution which will be appropriate m 
shallow water or where bottom interaction will be important. 
They further state that the FD and ODE methods will be 
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superior to SSFFT tecause an artificial bottom is not 
required for the benefit of the FFT. Explicit and implicit 
finite difference schemes are considered as well as an ordi- 
nary differential equation. The explicit finite difference 
(EFD) scheme is 
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where u n represents i(i(Rn,Zm), Rn is the nth range Dcint, and 
m 

Zm is the mth depth point. The implicit finite difference 
(IFD) method from Lee and Botseas (1982) is 
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The FD scheme is consistent, stable and convergent in the 
solution of the ?E (Lee and Papadakis, 1 979). 

The CDE method is 



d u = a u + b«, (u_ , , -2u_ +u_ , ) 

— m m m m + l m m - 1 * 

dR (AZ) 2 



where a and b are related by 



u_ = a(K , R , Z ) u + b(K ,R,Z)u„ 
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This technique is stable, consistent and convergent also 
which supports its use. One additional aspect of the ODE 
approach is the calculation of a variable range step. 

2. Compa ri son of PD and ODE w^th PE 

One result of the FD and ODE methods is greater 
flexibility in arbitrary bottom and surface boundary condi- 
tions than that available with the 5SFFT. Lee and Papadakis 
(1979) state t hat while the FD and ODE methods are both 
superior to the SSFFT, the best appears to be the ODE, fol- 
lowed by the EFD , which is limited by step size. Thus the FD 
and ODE methods appear to be significant improvements over 
the SSFFT. This result should extend to a WT based upon FD 
or ODE input. 

Wore recently Lee et al (1981) have re-examined the 
IFD and found it to be as equally efficient as the ODE. The 
general advantage cf the IFD over the SSFFT is that the 
problem is solved within the wafer column since the bottom 
effects can be modeled by straight line segments. This is 
similar to the ODE, but the IFD does not require the amount 
of storage necessary for the current ODE (Lee et al, 1981) . 
KcDaniei and Lee (198 2) have extended the IFD to the treat- 
ment of vertical density discontinuities. They found that 
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while without, this more realistic treatment of the water- 
bottom interface the IFD approximates the SSFFT and differs 
significantly from normal mode, the IFD corresponds 
extremely well with the normal mode calculations with the 
interface. In conclusion, the numerical techniques related 
to FD and ODE appear to be viable alternatives to the SSFFT. 
They will permit a more accurate representation of underwa- 
ter acoustic transmission, which will be directly relatable 
to the wt. 
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IV. ANALYSIS OF WT 
A. GENE HAL DESCRIPTION 

Twenty scenarios were selected to examine the general 
aspects of variations in several oceanic and geometric 
parameters. A set of input parameters and a sound speed 
profile were constant in each scenario while the source 
depth was varied. Five to ten source depths were used for 
each situation. The vertical and horizontal axes of the 
wavenumber spectra correspond to the spectral intensity and 
-he scaled wavenumber beta, respectively. The intensity 
axis was normalized to unity for ease in plotting. The null 
spacing was determined for each spectrum. These null spac- 
ing values were plotted versus source depths to yield the 
determination curve fcr each scenario. 

Table I lists the parametric specifications for each 
scenario. 3asically they can be divided info isospeed or 
sound speed profile cases. An isospeed profile of 4890 ff/s 
from the surface to the bottom was selected to evaluate the 
transmission geometry and the general consistency of experi- 
mental results, assuming a f ull y-absorbir.g bottom. The WT 



62 



is based on the arrivals of the direct and surface-reflected 
rays; thus, the isospeed cases will be the easiest to com- 
pare with theory. There will be no sound speed differences 
as a function of depth which could modify the arrival times. 
Additionally, several realistic ocean profiles were studied 
to allow a qualitative comparison with theory and explore 
spatial and temporal variability. It is recognized that the 
theoretical Lloyd's Mirror effect applies only to isospeed 
profiles. Results from real profiles were compared with this 
theory to see if they could approximate it. The spatial and 
tempera! scales of variability used in this research were 
approximately 200 miles and 20 days. 

More analyses were performed with the isospeed cases 
since this was an initial investigation. It still remains 
to consider ether simple cases, such as Linear negative and 
positive gradients, surface ducts and deep sound channels. 
Also realistic ocean studies that would involve greater tem- 
poral and spatial variations would be beneficial. Time was 
limited during this research, thus efforts were divided 
between model validation and oceanic sensitivity. It is 
believed, however, that the comparisons which will be 
presented are sufficient for a preliminary investigation 
into the sensitivity cf the WT. 
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TABLE I 

Scenario Specifications 



ISCSPEED PROFILE 



R un 
Set 


Rev Depth 
RD: ft 


Range 
R :nm 


Frag 
F : Hz 


Hater Depth 
Z: ft 


Number of 
SD Runs 


10 


300 


50 


100 


10000 


7 


11 


300 


50 


50 


10000 


6 


12 


300 


100 


50 


10000 


6 


13 


300 


50 


100 


5000 


6 


14 


300 


100 


100 


5000 


5 


15 


300 


50 


50 


5000 


6 


16 


300 


100 


50 


10000 


6 


17 


300 


100 


100 


10000 


5 


18 


800 


50 


50 


10000 


5 


19 


10000 


50 


50 


10000 


5 


R 1 


300 


25 


100 


10000 


5 


R 2 


300 


50 


100 


10000 


5 


R 3 


300 


75 


100 


10000 


5 



REALISTIC PROFILE 



R un 
Set 


Rev Depth 
R D : f 


Range 
R :nm 


Freg 
F : Hz 


Hater Depth 
Z: ft 


Src Depth 
SD (number) 


A0 


300 


100 


100 


9000 


8 


A 1 


1000 


100 


100 


9000 


5 


A 2 


9000 


100 


100 


9000 


5 


A3 


300 


100 


100 


9000 


5 


B0 


300 


100 


100 


12000 


3 


CO 


30 0 


100 


100 


9000 


8 


DO 


300 


100 


100 


12000 


8 
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The complex pressure input data were not modified so 
that the entire wavenumber spectrum would be observed at the 
receiver. These effects would need to be evaluated if the 
KT was to be compared to received signals. The surface and 
bottom were assumed flat. Boundary variations would have to 
be considered for WT use under realistic ocean conditions. 
A nominal beam width of thirty degrees was used and in 
almost all cases the range step was 0.1 mile. 

Frequencies of 50 and 100 Hz, water column depths of 
5000 and 10,000 feet and ranges of 50 and 100 nautical miles 
were the basic comparisons for most isospeed scenarios. A 
receiver dep^h of 300 feet and a fully absorbing bottom were 
utilized, unless otherwise noted. Additional comparisons 
were performed, such as range variations to 25 and 75 miles 
and receiver depth variations of 800 and 10,000 feet. A 
range of 100 nautical miles and a frequency of 100 Hz were 
used for the actual cases unless otherwise noted. The bottom 
depths were SOOO feet for profiles A and C and 12,000 feet 
for profiles B and D. 

The realistic profiles were obtained from the Acoustic 
Storm Transfer and Response Experiment (ASTREX) which was 
conducted by NPS in the northeast Pacifio Ocean in November 
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and December, 1980 (D unlap, 1 982} . The ASTREX study was con- 
ducted to evaluate the effect of oceanic response to atmos- 
pheric storms on acoustic propagation. The profiles were 
selected for this research to represent only spatial and 
temporal variability without respect to any experiment oper- 
ations. The profiles were obtained from expendable bathy- 
thermographs (AXBT) which were dropped by 0. S. Navy P-3 
aircraft. Fig. h shows the relative locations and times for 
these profiles. The digitized profiles were analyzed by the 
University cf Hawaii to remove spikes and false values. 
These profiles were combined with climatological salinity 
profiles at the Fleet Numerical Oceanography Center in Mon- 
terey. Table II gives the descriptive data for each 

profile. The bottom loss curves which were used are similar 
to those found in Urick (1979) and were arbitrarily devel- 
oped to indicate higher loss for profiles 3 and D. This is 
related to bottom composition; bottom types for profiles A/C 
and 3/D are calcareous sand and clay, respectively 
(NA VCCEA NO , 1 978) . Thus higher loss is expected for B/D. 

A general feature of the spectra for almost all rest 
cases is a distinct (J-shape wirh peak values at the minimum 
and maximum wavenumbers. Figs. 5 through 10 show this phe- 
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cure a. Space and time relationships of actual profiles. 
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TABLE II 

Location Information for Realistic Profiles 



Profile 


Lat 


Long 


Da te 


Time (Z) 


MLD (m) 


A 


41-24N 


127-23W 


16 NOV 


1928:04 


32 


B 


45-47N 


133-05W 


16 NOV 


2121:25 


61 


C 


41-23N 


128- 22W 


6 DEC 


1843: 37 


58 


D 


4 5-46 N 


133-09W 


5 DEC 


2033:45 
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r.omenon for isospeed cases and 
realise! c sound speed profile, 
figures are source depth (SD) , 
water column aep-.h (Z) in feet; 
and frequency (F) in Hz. This 
source depth determination curv 
later. The righ^ and left inten 
spond to beam e!evea*ion angles 
-ively, which could be consider 
ar.d surface-reflected waves. A 
thus, there may be a relation 
This beam angle is also near the 
(1 977) considered to be appropri 
might be a result of algorithmic 
angles could be tested to deta 



Figs. 11 through 16 for a 
The abbreviations on these 
receiver depth (RD) , and 
range (R) in nautical miles; 

key will also apply to the 
es which will be presented 
sity maxima appear to corre- 
ct 0 and 30 deg, respec- 
ed to approximate the direct 
30 deg beam angle was used; 
ship with the left maximum. 

maximum angle which lappert 
ate. Then the l a ft maximum 
inaccuracy. Different beam 
rmine if this would explain 



63 



the unexpected tJ-shace curve. This apparently anomalous 
spectral pattern might be explained by a possible transform 
assumption of symmetry about zero and the inclusion of neg- 
ative transform values. It is also possible that the left 
wavenumber peak may be the result of aliasing of the high 
wavenumber values. DiNapcli (1971) suggested that aliasing 
could occur and thus he utilized only one-half of the range 
field. This has not been tested in the current research; 
however, examination of different range segments might 
resolve this apparent problem. A remote possibility is that 
the direct and sur face-reflected rays eaoh exhibit a Bessel 
function display. This might be depicted by a peak at the 
reference wavenumber for each ray with a set of interfering 
envelopes between the peaks. A mathematioal explanation for 
this possibility has not yet been explored. No spikes ar c 
seen in very shallow water as depicted in Fig. 17, but also 
no distinct null pattern is observed, thus near surface 

source depths have been neglected. 

These figures were also included to show a series of 
source depth runs and the resulting decrease in beta spacing 
as a function of source depth. The null spacing was manu- 

ally measured for each run as the distance between two adja- 
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NORMALIZED SPECTRAL INTENSITY 




. Iscsoeea 
200 ft. 



re 5 



Profile Wavenumber Spectrum for SD 



NORMALIZED SPECTRAL INTENSITY 

O.OOO. DOS 0.010 O.OIS 0.020 0.023 0.030 0.035 0.0*0 0.0*5 0.050 



RUN M4 




Figure 6 . 



IsosDeed Profile Wavenumber Specirum for SD of 
300 ft. 
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i NORMAL 1 ZED SPECTRAL INTENSITY 

O.OOCD.OK O.OIO 0.013 0.020 0.023 0.030 0.033 0.040 0.043 0.030 



RUN 145 




gure 7. lecsoeea Profile Wavenumber Soactrum for SD of 
500 ft. 
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NORMALIZED SPECTRAL INTENSITY 

O.OOO.OOS O.QIO 0.01$ 0.020 0.02$ 0.030 0.03$ 0.010 0.015 O.OSO 



RUN 146 




Figure 3. 



IscsDeea Profile Wavenumber Spectrum for SD of 
300 ft. 
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NORMALIZED SPECTRAL INTENSITY 

O.OOQD.OOS O.OJO O.OIS 0.020 0.02$ 0.030 0.03S 0.010 0.01$ 0.0$0 



RUN H7 




. Is c= Deed Profile Wavenumber SDsetrura for SD o 
1505 ‘ft. 



7 '4 



gure 9 



NORMALIZED SPECTRAL INTENSITY 

O.OO®. 005 0.010 O.OIS 0.020 0.02S 0.030 0.035 0.040 0.015 0.050 



RUN 148 




cure 10. 



IscsDeed Profile Wavenumber Specrrum for SD of 
3000 ft. 
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NORMALIZED SPECTRAL INTENSITY 



RUN R03 

S0:200FT RD:300FT FrlOOHZ Z:9000FT R:100NM 




gure 11. 



Healisxic Profile Wavenumber Spectrum for SD of 
200 ft. 



75 



NORMALIZED SPECTRAL INTENSITY 

0.2 0.3 0.1 O.S 0.6 0.7 O.B 0.9 



RUN flCK 

SD:300FT RD:300FT FslOOHZ Z:9000FT R:100NM 




0.050 0.015 0.010 0.03S 0.030 0.025 0.020 0.0)5 0.010 0.005 0.000 

SCALED WAVENUMBER BETA (1/FT) 



gure 1 2. 



Real isric Profils Wavenumber SDectrum for SD o 
300 ft. 
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NORMALIZED SPECTRAL INTENSITY 

0.2 0.3 0.1 0.6 0.6 0.7 0.6 0.9 



RUN R05 

SD:SOOFT R0:300FT FilOOHZ Z:9000FT R:100NM 




gure 13. 



Realisxic Profile Wavenumber Spectrum for SD o 

550 ft. 
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NORMALIZED SPECTRAL INTENSITY 

0.2 0.3 0.1 0.5 0.6 0.7 0 



RUN R06 

SD:800FT RD:300FT F:100HZ Z:9000FT RslOONM 

O 




gure iu. Realistic Profile Wavenumber Soectrum for SD of 
800 ft. 
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NORMAL I ZED SPECTRAL INTENSITY 



SO: 15Q0FT RD:300FT 

3 



RUN fl07 

F : 100HZ Z:9000FT RtlOONM 



« 

Q 




i r ™r * ,- " r r i ■ — i — “ — i 1 

0.050 0.045 0.040 0.035 0.030 0.02S 0.020 0.015 0.010 0.005 0.000 



SCALED HAVENUM8ER BETA (l/FT) 



ire 15. ?£^ii§tic Profile 



Wavenumber Spectrum for SD of 
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NORMAL I ZEO 5PECTRAL INTENSITY 



RUN R08 




gune 16. Realistic Profile wavenumber Soecrrum for SD o 
3C00 ft. 
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NORMALIZED SPECTRAL INTENSITY 



RUN 141 

SD:50FT RD:300FT FrlOOHZ Z:5000FT RilOONM 




gure I"?. 



Iscsceea Profile Wavenumber SDactrum for SD o 
50 ft. 
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cent intensity minima, such as shown ir. Figs. 5 through 16. 
The nulls selected for measurement were those in the middle 
of the spectrum, which led to a consistent technique for 
comparison and avoided the rapidly decreasing null spacing 
near the left intensity maxima. The distance was measured 
to the nearest 1/64 of an inch and equated to the respective 
null spacing. Thus, there is some inhere?.* measurement 
error. The exact null locations were difficult tc determine 
in the shallow isospeed cases because of gen* le spectrum 
curvature and uncertainty in null location. The actual 
prefiles were very complex and again the nulls were diffi- 
cult to locate. One value for each spectrum plot was used 
to generate a point on the determination curve, such as Fig. 
19 fer a series of rcis. The variation of the data points 
from the theoretical Lloyd's Mirror curve will be the normal 
distance between the point and the curve. Measurement inac- 
curacy will decrease as source depth increases. This will 
be significant because at shallow depth there is some varia- 
tion noted from sound speed differences. The nulls are 
well-defined at areat depths, but there is lit*le variation 
already, thus the measurement errors were net considered 
significant. It is possible that a numerical algorithm which 
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could quanti- atively identify the null locations could mini- 
mize this error. 

This expectation is inferred from the wavenumber plots 
by Lauer (1979) and his examination of range determination. 
He found that range determination resolution increased as 
the distance between source and receiver decreased. It is 
possible that, a similar relationship might exist for depth 
determination. Thus greater resolution may be realized at 
shallow source depth. Even at greater source depths the res- 
olution may be acceptable for certain applications where an 
approximate depth is acceptable. 

The prediction curves which result from the individual 
WT spectra will be discussed in following sections. Table 
III is a listing of the prediction curve comparisons which 
were conducted. In the figures only one symbol may appear 
at a source depth if both curves had the same null spacing 
value. 
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SOURCE DEPTH (FT) 

0.0 5G0.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 



Figure 



SOURCE DEPTH DETERMINATION CURVE 




BETA NULL SPACING U/FT) 



13. ExamDle of Source Depth Determination Curve. 
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TABLE I II 

Comparison Specifications 



ISOSPEED PROFILE 



Figure 


Run Sets 


Constant Parameters 


Varying Parameter 


19 


11,15 


R: 50, F :50, RD: 300 


Z : 1 0000 ,5000 ft 


20 


10,13 


R : 50, F: 100, RD :300 


Z : 10000, 5000 ft 


21 


12,16 


R: 100 , F : 50 , RD :300 


Z: 10000,5000 ft 


22 


17,14 


R: 100, F: 100, RD: 300 


Z : 10000, 5000 ft 


23 


15,16 


F: 50, Z : 5000, RD: 300 


R : 5 0 , 1 00 nm 


24 


11,12 


F: 50, Z : 10000, RD:300 


R : 5 0, 1 00 nm 


25 


R 1 # R2,R3 


?: 100, Z: 10000 , RD: 300 


R:25,50,75 nm 


26 


12,17 


R : 100 , Z : 10000 , RD : 300 


F: 1 00, 50 hz 


27 


10,11 


R: 50, Z: 10000, RD:300 


F: 1 00,50 hz 


23 


11,18,19 


R : 50, F:50, Z: 1 0000 


RD: 300,800, 10000 ft 


REALISTIC PROFILE 






Figure 


Rur. Sets 


Constant Parameters 


Varying Parameter 


29 


AO, BO 


R: 100 , F : 1 0 0, RD:300 


S pace: 200 nm 


30 


CO, DO 


R: 100, F: 100, RD: 300 


Space:200 nm 


31 


AO, CO 


R : 100 , F: 100, RD: 300 


Time:20 days 


32 


B C , DO 


R: 100 , F: 100, rD: 300 


rime:20 days 


33 


AO, A3 


F: 100, 16N07, RD: 300 


R : 1 00 , 50 nm 


34 


AO, A1, A2 


F: 100, 16N0V, R:10Q 


RD: 300 , 1000,9000 ft 


COM PCSIT 


z 






Figure 


Run Sets 


Constant Parameters 


Varying Parameter 



35 (All Isospeed and Actua 1 Profile Run Sets) 
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B. ISCSPEED CASES 

1 • Wa ter Co lump. Ce £t h Vari a tion 

Fig. 19 is the comparison for a frequency of 50 Hz 
and a range of 50 nm at water column depths of 5000 and 
10,000 feet. It was not expected that there would be any 
variation if the water cclumn depth was modified, since all 
bottom interacting rays should be fully absorbed. The 
curves showed excellent agreement at source depths of 800 
feet and greater, but there was some variation noted for 500 
feet and shallower. Some variation may be related to meas- 
urement error, as noted above. Maximum variation from 

~heory is observed for both depths in the middle depth 
range, approximately 500 feet, which is also the region of 
greatest curvature in the source depth determination curve. 

Fig. 20 is the same except the frequency has been 
increased to 100 Hz. There is better agreement from 200 to 
3000 feet; however, the significance of the frequency varia- 
tion has not been evaluated. Fig. 21 is for a frequency of 
50 Hz and a range of 100 nm, and Fig. 22 is for 100 Hz and 
100 nm. Some variations are noted; however, they may not be 
significant, which would agree with the supposition that 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED WATER DEPTH VARIATION 
R:50NM F:50HZ R0:300FT 



o 




i- 



o 




0.0000 



LEGEND 

• - LL0Y0S MIRROR 

• - 11: Z-10000FT 
o-I5: Z-5000FT 




I ! I I I I 

0.0025 0.0050 0.0075 0.0)00 0.0125 0.0150 

BETA NULL SPACING (1/FT) 



0.0175 



Fiaure 19. IsosDe®a Frofile SD Dets rmination Curve: 
Variation - 1. 



00200 



*7 

U 
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water column depth differences will not be significant for 
an isospeed profile. 

2. Ra nge Variation 

Fig. 23 is the comparison for a frequency of 50 Hz 
and a water column depth of 5000 feet for 50 and 100 nm. 
There is essentially no difference between the curves. The 
experimental curves agree well from 800 to 3000 feet, but 
once again the shallow variations may result from inaccura- 
cies in evaluating beta from the spectrum plots. Fig. 24 is 
a similar comparison, but the water column depth has been 
increased to 10,000 feet. Fig. 25 increases the frequency 
to 100 Hz for ranges of 25, 50 and 75 miles. Again all 

curves show very good agreement with theory, which indicates 
that range variations within 100 nm of the source will prob- 
ably be insignificant. Ranges beyond 100 nm were not tested 
since Fitzgerald (1975) suggested that the SSFFT was not a 
good approximation for 100 HZ beyond 111 km. This could be 
tested later. 

3. Frequ e ncy Variation 

Fig. 26 is the comparison for a range of 50 nm and a 
water column depth of 5000 feet for 50 and 100 Hz. Both 
curves agree well with theory below 800 feet, but at 500 



39 



SOURCE DEPTH l FT) 

1000.0 1500.0 2000.0 2300.0 3000.0 3500.0 4000.0 



SOURCE DEPTH DETERM l NRT I ON CURVE 
ISOSPEED WATER DEPTH VARIATION 
R:50NM FslOOHZ RD:300FT 




LEGEND 

• - LLOYDS MIRROR 
■ - IOs Z-10000FT 
a - 13: Z-5000FT 



s- 



0.0000 




0.0025 0.0050 Q. 0075 0.0)00 0.0125 0.0150 

BETA NULL SPACING t 1/FTJ 



□ 



— r ■ — 1 1 

0.0175 0.0200 



gure 20. Isceoeed Frofile SD Determination Curve: Z 
Variation - 2. 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED WATER DEPTH VARIATION 




Figure 21. Iscspeed Profile 3D Determination Curve: Z 
Variation - 3. 
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SOURCE DEPTH DETERMINATION CURVE 
. ISOSPEED WATER DEPTH VARIATION 
. R: 100NM F : 100HZ RD:300FT 




gure 22. Iso^ne^d Profile SD Determination Curve: Z 
Variation - 4. 
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SOURCE OEPTH (FT) 

0.0 500.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 



SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RANGE VARIATION 
F:50HZ Z:5000FT RD:300FT 




LEGEND 

• - LLOYDS MIRROR 
■ - 15: R-50NM 
o - 16: R-100NM 



0.0000 




— i i 1 1 1 1 — 

0.002S 0.0050 0.0075 0.0100 0.012S 0.0150 

BETR NULL SPACING t 1/FT) 



— i 1 

0.0175 0.0200 



gure 23. 



Isosoeed Profile SO Determination Curve: R 
Variation - 1 . 
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SOURCE OEPTH (FT) 

0.0 SOO.O 1000.0 ISOO.O 2000.0 2500.0 3000.0 3500.0 4000.0 



SOURCE DEPTH DETERM I NRT I ON CURVE 
ISOSPEED RANGE VARIATION 
F :50HZ Z: 10000FT RD:300FT 



LEGEND 

• - LLOYDS MIRROR 
■ - 1 1 : R-50NM 
a - 12: R-100NM 




‘ » 

0.0000 0.0025 



~CUQ20O 



0.0050 0.0075 0.0100 0.0125 

BETA NULL SPACING ll/FT) 



0.0150 



0.0175 



gure 24. Isc^Deed Profile S D Determination Curve: R 
Variation - 2. 
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I SOURCE DEPTH l FT) 

0.0 500.0 1000.0 1S00.0 2000.0 2500.0 3000.0 3500.0 1000.0 



SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RANGE VARIATION 
F:100HZ Z :10000FT RD:300FT 




LEGEND 

• - LLOYDS MIRROR 

• - Rl: R-25NM 
a - R2: R-50NM 

• - R3: R-75NM 




1~~ i i i i i i 1 1 

0.0000 0.0Q2S 0.0050 0.0075 0.0100 0.0125 0.0150 0.0175 0.0200 

BETR NULL SPACING (1/FT) 



gure 25. 



Isospeed Profile SD Determination Curve: R 
Variation - 3. 
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feet the higher frequency appears to agree better. Fig. 27 
extends the water column depth to 10,000 feet and there is 
apparently better agreement with the theoretical curve. 
Overall there appears to be better agreement at 100 Hz than 
at 50 Hz, but this result is not certain with the limited 
number of cases run. 

4 . Recei ver Depth Variation 

Fig. 28 compares receiver depths of 300, 800 and 

10,000 feet for a range of 50 nm, water column depth of 
10,000 feet, and frequency of 50 Hz. Again all curves agree 
reasonably well. There is perhaps slightly better agreement 
for the shallow receiver, but this supposition cannot be 
verified from these limited data. 

5. Analysis of Isospeed Cases 

There do not appear to be any significant differ- 
ences for wazer column depth, range, frequency or receiver 
depth over the limits of variations considered. This was 

expected since the theory is supposed zo be only a function 
of source depth for an iscspeed profile. Thus “he WT would 
be expected to be equally effective for various transmis- 
sion geometry patterns and frequencies, if the sound speed 
profile is constant and the bottom fully absorptive. There 
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SOURCE DEPTH (FT) 

I 000.0 1SOO.O 2U00.Q 25£ 



source: depth determination curve 

ISOSPEED FREOUENCT VARIATION 
RtlOONM Z: 10000FT R0:300FT 
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aure 26. Isc$oeed Profile SD Ders rminetion Curve: F 
Verier ion -1. 
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SOURCE DEPTH (FT) 

1000.0 1300.0 2000.0 2300.0 3000.0 3300.0 1000.0 






SOURCE OEPTH DETERM I NRT I ON CURVE 
ISOSPEED FREQUENCY VAR I FIT I ON 
R:50NM Z:10000FT RD:300FT 
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Figure 27. 



Isc$oeed Profile SD Determination Curve: F 
Variation -2. 
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SOURCE OEPTH DETERMINATION CURVE 
ISOSPEED RCV DEPTH VARIATION 
R:50NM F:50HZ Z:10000FTi 




gure 23. 



Iscspeed Frofile SD Determination Curve: SD 
Variation . 
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may be slight effects for receiver depth and frequency, but 
these have not been substantiated. 

C. REALISTIC PROFILE CASES 
1 . S pat ial Variation 

Fig. 29 is the comparison of spatial variability for 
the curves which resulted from the spectra of profiles A and 
B. (The sound speed profiles corresponding to the pre- 
diction curves will be included for each comparison.) These 
were both sampled or. November 16 , 1980, and represent a spa- 
tial separation of about 200 nm. Both profiles result in a 
slight bias below the theoretical Lloyd’s Mirror curve. 
There is more variation along each curve, which is to be 
expected because the profiles are not isospeed. Also, par- 
tially absorbing bottoms are now used as compared to the 
fully absorbing bottoms used with the isospeed cases. The 
profile and bottom effects have not yet been separated to 
determine which contributes more to the departure from 
theory. This would require an extensive set of intercom pari- 
sons between these two parameters. (The present intent is 
to show that ocean variability is important.) It is diffi- 
cult to differentiate between the two curves, but it appears 
that there is a slight trend in 3 toward smaller null 
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spacings. Sound speed profile B is slowar than A over the 
entire depth. A slower sound speed would appear to indicate 
that a ray of lesser angle would be required to provide the 
second arrival for the interference pattern. That would cor- 
respond to a greater vertical wavenumber and decreased beta 
null spacing. 

Fig. 30 compares profiles C and D, which were sam- 
pled at the same locations as profiles A and 5, respec- 
tively, but on December 6, 1980. A similar pattern is noted 
with D spacing less than C, which may again correspond to 
slower sound speeds for profile D. 

2. Tempo ral Variation 

Fig. 31 compares profiles A and C which represent 
the same location, tut a separation of 20 days. There 
appears to be a curve crossing at about 800 feet, such that 
spacing is less at A for shallow depths ana greater for A 
below 300 feet. The A sound spped profile is slower than the 
C profile down to 800 feet and the pattern is reversed from 
800 to 2000 feet. Slower speeds appear to be related to 
shorter null spacings, which appears consistent with this 
example. The variations between the curves, however, are 
quite small. 
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SOURCE DEPTH DETERM I NRT I ON CURVE 
REALISTIC SPATIAL VARIATION 
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gure 29. 



Spatial Variation (A ,3 ) : SD Determination and 

AXBT Curves. 
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Fiaure 30. Soatial Variation (C ,D) : SD Da termination and 

AXBT Curves. 
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aure 3 1. Temoorai Variation (A,C) : SD Determination and 
A a 3T Curves. 
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Fig. 32 presents a similar comparison for profiles B 
and D. The speed profiles show almost no differences, except 
that D is slightly slower. This result appears as a 
slightly sharper deep sound channel which will mean slightly 
slower velocity for the direct path. That effect may be 
reflected as a slightly greater null spacing as shown. 

3. R ang e Va riat ion 

Results for profile A wer c compared for ranges of 50 
and 100 nm in Fig. 33. No significant variations were 
noted, which is similar to the results for - t he isospeed 
cases. 



4 . Re cei ver Depth Var ia tion 

Results for profile A were also compared for 

receiver depths of 300, 1000 and 9000 feet in Fig. 34. 

There may be slightly better agreement for the shallow 
receiver, but this is not certain. 

D. COMPOSITE COMPARISON 

Fig. 35 is a composite of all isospeed and actual 
profile data sets. Isospeed data show slightly better 
agreement with theory; however realistic profiles are a rea- 
sonable approximation to the theoretical curve. Realistic 
profile variations from the isospeed curve appear to be 
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related tc the sound speed profile and different bottom 
boundary conditions. The general relationships of the data 
indicate than below about 800 feet no significant variations 
are noted. Specific curves will apparently provide better 
predictions at shallower depths. It may be sufficient for 
certain applications no provide climatological or theoreti- 
cal prediction curves. These curves could be enhanced or 
modified by updated curves that would reflect the current 
sound speed profile. 
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gure 32. Temporal Variation (3,D) : SD Determination and 
AX3T Curves. 
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gure 33. Range (R) Variation for Profile A. 
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gure 35 



V 



CON cm 51 ON s 



The Wavenumber Technique ( WT) , as proposed by Lauer 
(1979) has been investigated for source depth determination. 
Several idealized and actual scenarios of oceanic and acous- 
tic parametric variations were tested. There is excellent 
agreement between the ideal test cases and the theoretical 
Lloyd's Mirror effect. The realistic cases showed signifi- 
cant similarity as well, but with some bias toward smaller 
wavenumber spacing. This is qualitatively attributed to 
sound speed profile variations in the water and sediment. 

There appears to be a response of the WT to variations 
in input, which is beneficial in allowing evaluation of the 
relative importance of individual parameters. Changes in 
range and water column depth do not indicate any significant 
variation. Some differences are possibly suggested for 

receiver depth and frequency, but the results have not been 

• 

fully analyzed. The WT appears to respond well to changes 
in the environmental parameters, sucn as sound speed 
profile. This is seen from the average bias of the actual 
cases from the isospeed cases. The bias does appear to be 
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relatively consistent and shows the importance of oceanic 
variations in both space and in time. Thus, it appears that, 
if information is available for the bottom structure, and if 
a realistic sound speed profile can be generated, then a 
very specific family cf curves could be generated for vari- 
ous frequency and receiver depth combinations. This simplic- 
ity indicates possible operational applicability and 
flexibility. 

These are the first known source depth determination 
curves for the WT and thus there is no comparison available. 
It may be that there are some biases in these curves which 
would be discovered when comparative runs are made using 
other inputs or propagation models. Comparison with observed 
signals might also reveal biases that are related to the WT 
or to the propagation model. An explanation has not been 
determined for the 0- shape nature of the wavenumber spectrum 
plots and thus slightly different curves might result if 
this is found to bias the null spacing. Also there is a 
bias associated with the method of manually estimating the 
null spacing. This might be resolved with a computer proce- 
dure to locate the nulls. The experimental curves actually 
agree well with theory under the assumption of a flat, 
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tot ally- reflecting surface and a fully-absorbing bottom. 
The premise that the theory might be approximately accurate 
in realistic cases appears valid. This overall generally 
smooth fit supports the validity of the wr. 

This has been an initial investigation of the WT, and 
thus only a broadbrush study was possible. These prelimi- 
nary results indicate that additional analysis might be ben- 
eficial. The response of the WT to input from different 
propagation models would determine if results from different 
models could be compared. Detailed investigation of the spe- 
cific effects of various parameters would lead to a quanti- 
tative analysis of the sensitivity of the WT. An explanation 
for the amomalous U-shaped spectrum and che development of 
an algorithm which can accurately determine null locations 
would be necessary before the WT could be used practically. 
An extremely important investigation would be the comparison 
of analyzed actual signals to predictions from the WT to 
determine applicability, accuracy, consistency and any 
appropriate biases. 

The intent of this research was to determine the valid- 
ity of the WT by development of source depth determination 
curves and to evaluate the oceanic and acoustic 
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ssnsitivit ies of these curves. The results obtained indicate 
that in general the WT appears to be valid wi^h respect to 
-he Lloyd’s Mirror effect. The parametric variations which 
were investigated show that oceanic variability is extremely 
important, but that the WT is flexible for other parameters, 
such as transmission geometry. The WT appears to be viable 
for operational applications. The additional investigations 
suggested above would determine the validity of this suppo- 
sition. It is conceivable that the WT could be developed 



into an operational prediction product. 
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